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A continuous wave atom laser formed by the outcoupling of atoms from a trapped Bose-Einstein 
condensate (BEC) potentially has a range of metrological applications. However, in order for the 
device to be truly continuous, a mechanism to replenish the atoms in the BEC is required. Here we 
calculate the temporal coherence properties of a continuously pumped atom laser beam outcoupled 
from a trapped Bose-Einstein condensate which is replenished from a reservoir at finite temperature. 

We find that the thermal fluctuations of the condensate can significantly decrease the temporal 
coherence of the output beam due to atomic interactions between the trapped BEC and the beam, 
and this can impact the metrological usefulness of the device. We demonstrate that a Raman 
outcoupling scheme imparting a sufficient momentum kick to the atom laser beam can lead to a 
significantly reduced linewidth. 

PACS numbers: 03.75.Pp, 03.75.Kk, 03.75.Dg 


I. INTRODUCTION 

A rudimentary atom laser can be realised by outcou¬ 
pling atoms from a trapped Bose-Einstein condensate. 
This has previously been achieved using either radio¬ 
frequency (RF) m 01 Raman transitions @,11 to trans¬ 
fer the BEC atoms to untrapped states which then prop¬ 
agate freely. The majority of experiments performed to 
date have outcoupled atoms from a BEC that is not re¬ 
plenished, leading to an overall decrease of beam flux 
as time proceeds. This places an upper limit on the in¬ 
tegration time for any measurement [7j. A truly con¬ 
tinuous wave (CW) atom laser will not only extend the 
integration time of any measurement, but has other po¬ 
tential advantages such as a reduction in the linewidth of 
the beam due to gain narrowing @], However, this will 
require a method of replenishing the number of atoms 
in the BEC so that the system achieves a steady state 
output. Continuous atom lasers have the potential for 
applications in metrology, interferometry and precision 
measurement 0113 , and the improved first-order coher¬ 
ence properties of BECs EMU may provide a number of 
advantages over thermal atom sources [T~il currently used 
in, for example, Sagnac interferometry for the detection 
of rotations [§, [ill |l6[ • 

The first experiment that demonstrated a continu¬ 
ous source of Bose-condensed atoms was performed by 
Chikkatur et al. G2- This was achieved by periodically 
replenishing a BEC held in an optical dipole trap with a 
new condensate transferred using optical tweezers. While 
the team was able to demonstrate a BEC of more than 
10 6 atoms at all times, this approach has the limitation 
that the coherence of the BEC is disturbed every time 
the condensates are combined due to their independent 
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phases. An alternate approach was demonstrated by La- 
haye et al. [?J who performed evaporative cooling on a 
magnetically guided continuous atomic beam. They ob¬ 
served a significant increase in phase-space density, but 
did not reach Bose-Einstein condensation. More recently 
Robins et al. have demonstrated a pumped atom laser 
whereby atoms from one condensate acting as a source 
undergo a stimulated transition into a second conden¬ 
sate from which an atom laser beam was extracted using 
RF outcoupling [13 . Recently a continuous strontium 
BEC formed by laser cooling has been demonstrated m. 
which is potentially an excellent source for an atom laser. 
While these experiments have demonstrated the separate 
elements necessary for a continuous atom laser, they have 
yet to be combined in a single experiment. 

Many theoretical proposals for continuously-pumped 
operation have been made, using either evaporative cool¬ 
ing [201423 or spontaneous emission [23i425j |. A number 
of studies have already considered the properties of a 
CW atom laser at zero temperature [2p| 1.361. However, 
one of the most likely experimental routes to a true CW 
atom laser will involve replenishing the BEC via cool¬ 
ing from a reservoir containing thermal atoms so that 
the condensate mode is maintained in a steady state 
via Bose-stimulated collisions, similar to the experiment 
of Stellmer et al. 0. This will avoid the deleterious 
linewidth and heating effects associated with the merg¬ 
ing of initially separated condensates with independent 
phases E2- However, this will also mean that a CW 
atom laser will necessarily operate at finite temperature, 
and hence the effect of thermal fluctuations on the atom 
laser coherence should be accounted for. An initial study 
of such issues has been performed by Proukakis @3] us¬ 
ing a classical field model for the BEC replenished by 
an undepleted thermal reservoir. Here we build on this 
work by implementing a conceptually similar model but 
with specific improvements. In particular, we implement 
an appropriate energy cutoff in the harmonic oscillator 



2 


basis for the classical field, and propagate the atom laser 
beam in a plane wave basis well beyond the influence of 
the trapped BEC, as illustrated in Fig. [Q We focus on 
understanding the linewidth of the atom laser beam as 
a function of both the temperature and the momentum 
kick of the Raman outcoupling. 


II. MODEL OF A CONTINUOUSLY PUMPED 
ATOM LASER 


We model a trapped BEC at finite temperature us¬ 
ing the simple growth theory of the Stochastic Projected 
Gross-Pitaevskii equation (SPGPE) formalism [38M40t | . 
in which the energy damping of the scattering terms 
are neglected. This has previously been successfully ap¬ 
plied to model the formation of Bose-Einstein conden¬ 
sates from evaporative cooling [4lJ . Briefly, the second 
quantised Bose-field operator is divided at a cutoff en¬ 
ergy E cut into a high occupation (Nk 1) coherent 
region, in which interaction effects and thermal fluctu¬ 
ations are significant, and a low occupation incoherent 
region that can be treated as a thermal gas. The dy¬ 
namics of the entire coherent region can then be treated 
approximately as a classical field. For the problem to 
be tractable we treat the incoherent region as a unde¬ 
pleted reservoir at fixed temperature T and chemical po¬ 
tential g. We expect this is a reasonable approxima¬ 
tion in steady state if we imagine that we can replenish 
thermal atoms at a constant rate while undergoing con¬ 
tinuous cooling, e.g. |l9|, IHJ . The stationary incoherent 
region thus continuously replenishes the coherent region 
through Bose-stimulated collisions. The atom laser beam 
is formed by implementing Raman outcoupling from the 
coherent region, which depletes not only the condensate 
but also the non-condensed atoms, resulting in a ther¬ 
mally broadened linewidth for the atom laser beam. 

We consider a magnetically trapped 8 'Rb conden¬ 
sate containing approximately 10 5 atoms. We outcou- 
ple atoms into an atomic waveguide, as demonstrated 
by Guerin et al. [42j | , by transferring them to a magneti¬ 
cally insensitive state via a Raman transition which gives 
the outcoupled atoms a momentum kick in the positive 
x direction. The waveguide provides confinement in the 
radial direction. We choose experimentally reasonable 
trapping frequencies of {uj x ,uj±) = 2ttx (10, 200) Hz. For 
relatively tight transverse confinement, we neglect the 
dynamics in the radial direction and simulate the system 
in one dimension. 

The equations of motion for the trapped field i/j(x,t) 


and outcoupled atoms 4>(x,t) are 


dip(x) = --V (L +g 12 \(/}(x)\ 2 \)ip(x) - hfle lk ° x cj){x) 

n L 


dt 


+ V 


d(f>(x ) = — - 
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with 


n 2 d 2 1 

L = ~ 2m dx 2 + 2 mU}xX +9ll ^ x ^ ’ ( 3 ) 


where m is the atomic mass, and g and T are the chemi¬ 
cal potential and temperature of the reservoir. The recoil 
momentum and Rabi frequency of the Raman transition 
are hko and Q respectively, with a two-photon detuning of 
S. The operator V projects tf>(x) onto the truncated basis 
of harmonic oscillator eigenstates with energy e < E cut . 

The effective ID interaction strengths are gij = 
where crj_ is determined by integrating out the transverse 
dimension using a Gaussian variational anstaz, and a^- is 
the s-wave scattering length for collisions between states i 
and j. In this work we choose an = ai 2 = a 2 2 = 5.29 nm, 
which is approximately true for the IE = 1 ,nif = — 1 ) 
and l-F = 2,mj? = 0) states of 87 Rb [43j. The noise 
is complex Gaussian and has a non-vanishing correlator 
(dW*(x,t)dW^(x',t)) = 2"/S(x — x')dt. The quantity 7 
is the rate constant for the exchange of particles between 
the thermal reservoir and the classical field, for which we 
use the estimate 7 = 12 maf 1 fcsT/ 7 rS 3 [i3| . The steady 
state of the atom laser is independent of this quantity 
inside the weak outcoupling regime such that the con¬ 
densate is not significantly depleted. We chose H = 1 
Hz for all our simulations in order to stay in the weak 
outcoupling regime [45|. 

We simulate the classical field ik(x) us ing the harmonic 
oscillator basis as described in Ref. [3a Sg , and present 
results for a cutoff of E cut = 150 Huj x . We have increased 
the cutoff to E cut = 300^07 in other simulations and 
found no significant difference in our results. The beam 
4>{x) is simulated on a rectangular grid, and we make use 
of numerical grid interpolation for the coupling terms. 
Classical field atoms are transferred to the beam with a 
momentum hko , and we implement an imaginary absorb¬ 
ing potential at the edge of the atom laser beam grid. 

We initialize our simulations with a T = 0 condensate 
found from the solution of the time-independent Gross- 
Pitaevskii equation with the outcoupling turned off. The 
simulations are then run until we reach a steady state for 
the beam and condensate. We calculate observables by 
averaging 500 time samples of 64-95 independent trajec¬ 
tories, which gives converged results with relatively small 
statistical noise. Fig. |T] shows examples of the instanta¬ 
neous density of the trapped Bose gas and the atom laser 
beam once the system has reached steady state at (a) 
T = 0, and (b) T = 144 nK. 
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FIG. 1. (Color online) Typical classical field and atom laser 
beam profiles in steady state, (a) T = 0, (b) T = 144 nK. The 
instantaneous classical field density is shown in blue, whereas 
the instantaneous atom laser beam density is the black-dashed 
curve. The red dot-dash curve shows the time-averaged atom 
laser beam density for comparison. The atom laser beam 
density has been scaled by a factor of 2 x 10 3 so that is visible 
on the axes as the trapped classical field. Insets: Detailed 
view of the classical field region. In (a), the instantaneous 
atom laser beam density and time-average atom laser beam 
density lie on top of one another. For all plots: g = 102 
hco x , gu_/hu x xo = 0.0166 where xo = (h/miJx) 1 ^ 2 , D = 1 Hz, 
fiy = 0.0028, k 0 = 2.3 x 10 7 m _1 . . 


III. TEMPORAL COHERENCE OF THE 
OUTPUT BEAM 


The phase coherence of the atom laser beam can be 
quantified by the normalised equal position first-order 
temporal correlation function [47| 




(4>* (x, t)cj>(x, t + At)) 

V(Wx,t)| 2 )Mx 1 t + Ai)|V 


( 4 ) 
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Figure [2] shows |</ 1 '(a;, At)| for two different tempera¬ 
tures. Unsurprising, the coherence time decreases mono- 
tonically with increasing temperature for fixed /r. In this 
parameter regime we find that the atom laser coherence 
directly reflects the first-order temporal coherence of the 
classical field ip(x) at trap centre, which is plotted for 
comparison in Fig. [2] 

To gain physical insight into the decrease of phase co¬ 
herence in the condensate at finite temperature, we can 
approximate the evolution of the classical field by utilis¬ 
ing the Bogoliubov approximation [48'j. Expressing the 
classical field as 


ip(x,t) = e~ l ^ h 



+ E x)e ie ^ R bf) 

j> o 


( 5 ) 


where No is the condensate number, ipo (a?) is the conden¬ 
sate wave function, and Uj(x) and v 3 (x) are the ID Bo¬ 
goliubov quasiparticlc modes with energies ej measured 
relative to the chemical potential [49], An approximate 
form for the Bogoliubov modes in the ID Thomas-Fermi 
limit can be written as [50 ) 


/f(i) 




P 0 {x), 


( 6 ) 


where the spatial dependence has been normalised to 
the Thomas-Fermi radius of the condensate x = x/ Rtf 
where I?tf = (2/Vtow-c) 1 / 2 , and Pj{x) are Legendre 
polynomials. The functions ff{x) are related to the Bo¬ 
goliubov mode functions u{x) and v{x) via the relation 
ff(x) = u(x) ± v(x). The Bogoliubov modes have ener¬ 
gies €j = hx xy /j(j + l)/2 for positive integers j [50] in 
the one-dimensional Thomas-Fermi limit. 

We can then construct an estimate of the first-order 
correlation function g^(x,At) oc (ip*(x, t)il>(x, t + At)) 
by substituting the Bogoliubov expansion (0 and com¬ 
puting the sum numerically. We examine the temporal 
coherence at the centre of the condensate where outcou- 
pling is resonant. We assume each Bogoliubov mode is 
thermally occupied with quasiparticles according to the 
equipartition distribution 


N 3 EE (b*b 3 ) = - 2 -, (7) 

e 3 

while the ground state has an occupation equal to the 
condensate number Nq calculated from the one-body 
density matrix using the Penrose-Onsager criterion for 
condensation [Hj]. Assuming Nj » 1, we have 


(tp*(0,i)ip(0,t + At)) 


N 0 \vom 2 +J2 — 

j> o 
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/ V + e J 
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FIG. 2. (Color online) First-order temporal coherence of the 
atom laser beam and condensate quantified by At) for 

temperatures (a) T = 48 nK with No = 1.1 x 10 5 and (b) 
T = 384 nK with No = 0.72 x 10 s . For both plots — solid 
red line is the coherence of the condensate at x = 0; dashed 
blue line: coherence of the atom laser beam for x just outside 
the trapped gas; dash-dotted black line: estimate given by 
Bogoliubov approximation of Eq. lf8]l. All other parameters 
are the same as for Fig. [T] 


The estimate given by Eq. (0 is also shown in Fig. [2j 
where the correlation function has been normalised to 
compare with the SPGPE results. The estimate repli¬ 
cates the early decay of the numerical results, indicating 
that the mechanism for the loss of coherence observed in 
the condensate is the deplrasing of thermal fluctuations. 
The revivals of coherence present in the Bogoliubov esti¬ 
mate are not present in the SPGPE result, as the noise 
term of the SPGPE causes a random drift in the phase 
of each of the modes and preventing the rephasing. 


IV. INFLUENCE OF THE MOMENTUM KICK 
ON THE LINEWIDTH 

We now investigate the effect of the magnitude of the 
momentum kick on the linewidth of the atom laser at 
finite temperature. At zero temperature the momen¬ 
tum transfer to the atom laser beam of hko would have 
no bearing on the temporal coherence and hence the 
linewidth of the atom laser. At finite temperature we 


find that the temporal coherence of the trapped BEC is 
also independent of the value of fco- However, we find 
that this is not the case for the atom laser beam. In¬ 
stead, below a critical value of we find that for a 
given temperature the temporal coherence of the beam 
decays significantly faster than that of the trapped BEC. 

Defining the power spectrum of a field xi x o>t) as 


P x (xo , w ) 


dte lult x(xo ,t) 


(9) 


we define the spectral linewidth of the field as 


Aw x (i 0 ) = 


J dojP x (xo, w)w 2 — (/ dojP x (xo,oj)uj 


1 1/2 


( 10 ) 


This quantity is plotted for both the centre of the con¬ 
densate (blue asterisks) and the middle of the atom laser 
beam (red diamonds) in Fig. [3] for a temperature of 
T = 384 nK. It is clear that the atom laser beam ex¬ 
hibits a dramatic increase in linewidth for momentum 
kicks below ko = 6 x 10 6 m -1 . 

We have found that this excess linewidth is due to the 
thermal density fluctuations of the classical field that the 
atom laser beam experiences as it propagates through the 
trapped BEC. If the outcoupling momentum hk is suffi¬ 
ciently large, the density of the trapped classical held 
is essentially frozen on the timescale required for parti¬ 
cles to leave the region of the BEC, so the linewidth of 
the beam is determined by the spectral linewidth of the 
condensate at the spatial position where the outcoupling 
is resonant. However, for smaller momenta, and there¬ 
fore longer exit times, the outcoupled atoms experience 
significant spatio-temporal fluctuations in the effective 
spatial potential due to the classical Held, resulting in 
greater dispersion in velocity of the outcoupled atoms 
and a larger linewidth. A threshold momentum separat¬ 
ing the different linewidth behaviour may be determined 
via a classical argument. Particles travelling significantly 
faster than the thermal density fluctuations can ‘outrun’ 
the fluctuations, and thus experience no velocity broad¬ 
ening. Taking the sound speed c s = (gn/V o |^(0)| 2 /m) 1 / 2 
as an estimate for the characteristic speed of density fluc¬ 
tuations, we note that the threshold momentum for the 
excess linewidth is approximately twice this value, con¬ 
sistent with our argument. 

We have implemented a simple numerical model to 
confirm this explanation. We have simulated the dynam¬ 
ics of an ensemble of classical atom laser beam particles 
placed at the classical outcoupling position with initial 
momentum hko- In a zero temperature system the effec¬ 
tive potential experienced by the outcoupled atoms due 
to interactions with the BEC will be static. Therefore, 
outcoupling the atom laser beam from the centre of the 
trap will result in a momentum increase of the output as 
the atoms “roll down the hill” so their output momentum 
is hk = y/(hko) 2 + 2 mji j35| , but there is no broadening 
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FIG. 3. (Color online) Linewidths of the condensate and the 
atom laser as a function of momentum kick for T = 384 nK. 
Blue asterisks: condensate linewidth at outcoupling position. 
Red diamonds: atom laser beam linewidth. Black dots: Clas¬ 
sical model of the linewidth of Newtonian particles exiting 
the system experiencing the fluctuating classical field of the 
trapped BEC as described in the text. All other parameters 
are the same as for Fig. [T] 


of the initial distribution. However, at finite tempera¬ 
ture, on exiting the trapped BEC, the momentum distri¬ 
bution of the beam will be broadened due to interactions 
with the fluctuating density of the classical field, and we 
quantify it by standard deviation of the output energy 
distribution. The resulting linewidth is plotted in Fig. [3] 
in comparison with the results from our SPGPE model, 
where we can see they are in broad agreement. 

As the excess linewidth of the atom laser beam is due 
to density fluctuations of the trapped BEC, the effect 
will disappear if there are no interactions between the 
trapped BEC and the atom laser beam. We have verified 
this by setting g 12 = 0 in our simulations (while appro¬ 
priately adjusting the detuning S = — g/h in 

order to maintain resonant outcoupling in the centre of 
the condensate.) In this case, the linewidth of the atom 
laser beam reverts to the linewidth of the condensate for 
all values of the Raman kick momentum. It may be possi¬ 
ble to exploit an inter-component Feshbach resonance to 
achieve this in an experiment, significantly reducing the 
linewidth. However, this may be technically challenging 
in practice. 

A potential solution to reduce the atom laser linewidth 
for a fixed temperature and momentum kick is to alter 
the detuning such that the resonant position is closer to 
the edge of the condensate and the exit time of the beam 
is shorter. We do observe a monotonic decrease in the 
linewidth as the outcoupling position is moved towards 


Rtf, but it remains larger than that of the condensate. 
Beyond Rtf , the atom laser linewidth approaches that of 
the thermal reservoir, as expected — only thermal atoms 
are outcoupled. As a result, the linewidth of the atom 
laser beam is at least that of a thermal state for all values 
of the two-photon Raman detuning. We conclude that 
altering the detuning may decrease the overall atom laser 
linewidth, but cannot reduce it to that of the condensate. 
Furthermore, any reduction in linewidth is offset by the 
reduction in flux (for a fixed f2) associated with moving 
the outcoupling position away from the centre of the trap. 

Another strategy to reduce the excess linewidth ob¬ 
served for a given temperature and momentum kick may 
be to remove the waveguide and allow outcoupled atoms 
to instead fall under gravity in the radial direction, po¬ 
tentially reducing the exit time of the atom laser beam. 
For the parameters used in this work, the gravitational 
acceleration results in an exit time equivalent to an initial 
momentum kick of 4.4 x 10® m -1 . This value is below the 
threshold momentum, and so the linewidth of the atom 
laser beam will still be excessive. Additionally, removing 
the waveguide will degrade the transverse pro file of the 
atom laser beam [42, |5l| , as was seen in [3, [521455} . 

V. CONCLUSIONS 

We have modelled the temporal coherence and 
linewidth of a continuous atom laser at finite temper¬ 
ature. We find that when outcoupling with a large mo¬ 
mentum kick, the first-order coherence of the atom laser 
beam mimics the first-order coherence of the trapped 
BEC which is determined by thermal fluctuations. How¬ 
ever, below a threshold momentum kick related to the 
condensate speed of sound, there is a dramatic increase 
in the spectral linewidth of the atom laser beam. We 
have identified the origin of this as density fluctuations 
in the trapped BEC while the atom laser beam exits the 
condensate, further increasing the energy spread of the 
output beam. Our results will influence the design and 
optimisation of continuous atom lasers at finite temper¬ 
ature in the future. 
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